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Abstract. 

We derive functional renormalization group schemes for Fermi systems which 
are based on the two-particle irreducible approach to the quantum many-body 
problem. In a first step, the cutoff is introduced in the non-interacting propagator 
as it is commonly done in functional renormalization group based on one-particle 
irreducible vertex functions. The most natural truncation of the resulting infinite 
hierarchy of flow equations is shown to be fully equivalent to self-consistent 
perturbation theory. An earlier suggested alternative truncation strategy is 
considered as well. In a second step, the cutoff is introduced in the two- 
particle interaction. Again two truncation procedures are investigated, one of 
which was derived before. In the latter, the mean-field solution of the many- 
body problem is considered as the starting point of the renormalization group 
flow. We compare the performance and the required numerical resources for 
solving the coupled flow equations for all the approximate schemes by applying 
them to the problem of the quantum anharmonic oscillator. In a functional 
integral representation, this model has a formal similarity to the quantum many- 
body problem. The perspectives for applying the derived two-particle irreducible 
functional renormalization group approaches to zero- and one-dimensional systems 
of correlated fermions are discussed. 
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1. Introduction 

Over the past years the functional renormalization group (fRG) was established as a 
versatile tool to study low-dimensional interacting Fermi systems mm- The variant 
almost exclusively used up to now is based on the generating functional for the one- 
particle irreducible (IPI) vertex functions. In this approach, one derives an exact 
infinite hierarchy of coupled differential flow equations for the IPI vertex functions 
(self-energy and effective n-particle interactions), where the derivative is taken with 
respect to a flow parameter. Prior to practical calculations for any fermionic many- 
body model, a truncation of this set of coupled flow equations is required. Most 
truncation schemes were guided by perturbation theory in the flowing two-particle 
interaction rendering IPI fRG a non-perturbative approximate method which is 
controlled for small to intermediate interactions. It is one of the assets of fRG based 
approximation schemes that they are unbiased in the sense that barely any a posteriori 
knowledge of the physics of the problem at hand must be used when setting them up. 
This e.g. allows for a faithful analysis of the interplay of distinct emergent quantum 
many-body phenomena. 

The IPI fRG was used to investigate the competing ordering tendencies (e.g. 
magnetism and unconventional superconductivity) in the two-dimensional Hubbard 
model and variants of the latter to study the spectral and transport properties 

of inhomogeneous one-dimensional lattice models (quantum wires) falling into the 
Luttinger liquid universality class [T]i "''^sll as of quantum impurity (quantum dot) 
models [T]. Gompared to other RG approaches, fRG has the distinct advantage that 
it can directly be applied to microscopic lattice models and not only to low-energy 
effective field theories. This way, high-energy features and complex crossover behavior 
are captured. 

For problems of interest in which plain perturbation theory in the two-particle 
interaction of characteristic amplitude U is divergent in the infrared limit (e.g. in 
two-dimensional models or in translationally invariant Luttinger liquids), the flow 
parameter must be introduced such that it has infrared regularizing properties; it 
acts as a low-energy cutoff [!]■ In. other systems, low-order perturbation theory 
is regular but, due to logarithmic terms with prefactor C/", n S N, restricted to 
parameter regimes of infinitesimal extend. Then, IPI fRG can be used to resum the 
logarithms generically leading to power-law behavior with [/-dependent exponents 
[1]. In both cases, characteristics are utilized which are commonly associated to RG 
methods, e.g. the successive inclusion of energy scales from high to low and the 
resummation of logarithmic terms. However, IPI fRG was also used for problems in 
which perturbation theory is regular and not plagued by logarithms. For such, fRG 
should be viewed as a “renormalization group enhanced perturbation theory” in which 
certain types of diagrams are resummed [UllllSlle]. It is then possible to employ other 
flow parameter schemes than the introduction of an infrared cutoff. One can e.g. think 
of successively turning on the strength of the two-particle interaction during the flow 
(see e.g. Ref. B)- 

Recent attempts to abandon the realm of small to intermediate interactions within 
fRG include to set up IPI schemes in which certain aspects of the interaction are 
already included in the initial point of the RG flow 0191 Ho]; e.g. the idea of taking 
dynamical mean-field theory containing all the local correlations as the starting point 
of the fRG flow was put forward [3] . In another attempt for a certain quantum impurity 
model, namely the single-impurity Anderson model, lots of a posteriori knowledge of 
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the physics was built into the fRG procedure m- In the latter case, however, one 
gives up the strength of fRG of being unbiased. The IPI fRG was recently extended 
to study quantum wires as well as quantum dots in non-equilibrium. The driven non¬ 
equilibrium steady state [niiiiiiiin] as well as the non-equilibrium time evolution 
[16j were investigated. 

Despite its successes, present day IPI fRG suffers from several shortcomings. 
For general two-dimensional models, it still constitutes a challenge to formulate an 
unbiased purely fermionic scheme which allows to extend the flow into phases with 
spontaneously broken symmetry [T]. Phase transitions are driven by composite fields 
and not the bare fermions. Grudely speaking, these composite degrees of freedom are 
not properly dealt with in the standard implementation of IPI fRG. As a consequence, 
the RG flow diverges at a certain cutoff scale and has to be stopped. Therefore, 
only ordering tendencies can be read off but one cannot compute observables of 
the cutoff-free problem in the symmetry broken state. In recently developed novel 
schemes the flow is initialized including a small selected symmetry breaking term. In 
the parts of the parameter space in which the corresponding instability dominates, 
the flow is regularized and the cutoff can be removed allowing to investigate the 
symmetry broken phase dzi nHi m uni HD- Another way out is to introduce the 
bosonic composite degrees of freedom during the RG flow by partial ’’bosonization” 
[21 1211 1231 [Ml IMl HSl HZ]: however, again to the price of giving up treating all 
instabilities on equal footing. 

For zero-dimensional quantum impurity models, the above mentioned “flow to 
strong coupling” is not an issue and the truncated set of IPI fRG flow equations 
can be integrated down to vanishing cutoff. From the resulting expressions for the 
vertex functions, approximate results for observables of interest such as the ground 
state energy or the single-particle spectral function of the cutoff-free model can be 
computed. However, case studies revealed that going to higher orders in the truncation 
does not necessarily lead to the hoped for improvement of the results beyond the limit 
of weak interactions (for which the results systematically improve). E.g. for the 
single-impurity Anderson model, the local spin susceptibility shows Kondo physics 
in first order truncation [^ but loses this property in second order HIS]. Similarly, 
for the interacting resonant level model, expected power-law behavior of observables 
confirmed in first order truncation is lost in second order |15j . 

A general drawback of the IPI fRG for practical computations is that it is 
computationally costly to systematically include frequency dependence. This restricts 
the access to dynamical observables. A frequency dependence of the self-energy is 
only generated if the dependence of the two-particle vertex on its three independent 
frequencies (for a time translational invariant model) is kept [1]. This leads to an 
excessive growth of the number of equations to be consideredjf] Therefore, it was 
only for zero-dimensional quantum dot models with a few interacting degrees of 
freedom H E] or toy models such as the quantum harmonic oscillator with quartic 
anharmonicity H that frequency dependence up to second order could be included 
completely. This e.g. prevented the study of bulk Luttinger liquid physics in one¬ 
dimensional models of interacting electrons. Schemes to include a frequency dependent 
self-energy for two-dimensional systems were put forward in Refs. I2H1IM11IS1- 

Finally, it remains a challenge to fulfill Ward identities in truncated IPI fRG 
approaches beyond the order in which the considered set of flow equations contains 

X For an approximate treatment of the frequency dependence, see Refs. IIEIITI]. 
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all diagrams of plain perturbation theory [301 Enin]. Therefore, the relation between 
approximate IPI fRG schemes and so-called conserving approximations |3H E3j 
remains elusive ||] 

To circumvent the first mentioned drawback of IPI fRG for problems of two- 
dimensional fermions, it was suggested in Refs. [33 |33| to set up fRG schemes 
based on the two-particle irreducible (2PI) approach to quantum many-body physics 
[32l [33] • This would allow to better incorporate the composite degrees of freedom. 
Such a procedure was earlier hinted at in Ref. m- However, recent results obtained 
by dynamical mean-field theory for models with phase transitions [381 139] reveal 
unforeseen divergencies in the 2PI two-particle vertex. These might lead to obstacles 
for a 2PI fRG treatment of two-dimensional models. 

We are primarily interested in interacting quantum dots and wires and here study 
2PI fRG approaches to investigate their potential to overcome the second and third 
problem of IPI fRG. Furthermore, using the same 2PI framework for fRG schemes as 
in standard conserving approximations promises insights on the relation to these. 

In a first step, we derive four pair-wise related but different 2PI based fRG 
schemes. Two of them were introduced earlier |35L 136] . For these, we suggest 
modihcations and provide insights which were not emphasized so far. Furthermore, we 
allude to the relations between all the approximate approaches. Instead of studying a 
fermionic many-body problem, we then apply these methods in a comparative study 
to the problem of the quantum harmonic oscillator with quartic anharmonicity. In a 
functional integral representation of the partition function, the anharmonic oscillator 
has the same structure as the one of a fermionic many-body problem (although with 
real helds instead of Grassmann fields). The anharmonicity g corresponds to the 
amplitude of the two-particle interaction in a many-body problem. Considering this 
toy model, we thus reduce the numerical effort (due to the reduced number of degrees 
of freedom and the more rapid decay of propagators) without sacrificing the general 
structure of the equations to be solved. We emphasize that the formal complexity of 
the problem we study is comparable to that of the single-impurity Anderson model 
[3] . The anharmonic oscillator was earlier used to illustrate the fRG procedure and to 
test newly developed fRG schemes [isnuniiiiisi- It has the additional advantage 
that numerically exact results for observables of interest can easily be obtained. By 
comparing exact results to our approximate ones for the dependence of the ground 
state energy and the position fluctuations on g as well as for the frequency 

dependence of the self-energy at fixed g, we obtain an impression of the approximation 
quality. The 2PI fRG results are in addition compared to the ones of IPI fRG. We 
furthermore discuss the advantages and drawbacks of the different schemes in their 
application to fermionic many-body problems. 

We briefly study another fRG scheme which is not based on a 2PI generating 
functional but on one which is two-particle point-irreducible (2PPI). The use of this 
was suggested [33J |33 for ifs relation to density functional theory and applied to the 
anharmonic oscillator in Ref. [35] . The results of the 2PPI approach to the anharmonic 
oscillator are compared to those of 2PI. 

We identify two 2PI schemes in which a frequency dependent self-energy is 
generated consistently without having to deal with functions of three independent 
frequencies. Compared to the standard frequency dependent IPI approach this 
considerably reduces the numerical effort. Nevertheless, the accuracy of these two 

§ For recent considerations on the relation between fRG and Schwinger-Dyson equations, see Ref. m- 
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2PI schemes for the anharmonic oscillator is satisfactory. For one of them, it is 
even slightly superior to IPI fRG. It remains to be seen if the same holds for 
fermionic many-body systems. We furthermore prove that the most natural truncation 
of the 2PI flow equations with the flow parameter introduced via the free single¬ 
particle propagator exactly gives the well known self-consistent perturbation theory 
(conserving approximation), e.g. in lowest order the self-consistent Hartree-Fock 
approximation. 

Our paper is organized as follows. In Sect. we give the functional integral 
representation of the partition function of interest. The generating functionals of 
the 2PI formalism are introduced in Sect. In Sect. we present a first 2PI fRG 
approach which results from including the flow parameter in the free propagator C 
(as is mostly done in IPI fRG); we coin it C-flow. The two approximation schemes 
resulting from C-flow as well as their application to the anharmonic oscillator are 
presented in subsections. As a sequel of Sect. we present technical details of the 
pair propagator and the Bethe-Salpeter equation in Sect. which are crucial for the 
second 2PI fRG approach, the so-called [/-flow. In this approach, the flow parameter is 
introduced via the amplitude U of the two-particle interaction as described in Sect.[^ 
Again, subsections are devoted to two resulting approximation schemes. In one, the 
non-interacting problem is considered as the starting point of the RG flow; in the other, 
the flow is started from the mean-field solution of the problem at hand. The results 
of these schemes for the anharmonic oscillator are presented in another subsection. 
In Sect. [Tj we give a brief account of the 2PPI fRG and discuss our results for the 
anharmonic oscillator. The results for the most successful approximate 2PI and 2PPI 
schemes are compared to the ones of IPI in a concluding Sect.[^ It also contains a brief 
discussion of the perspectives of applying 2PI fRG to fermionic quantum many-body 
problems. In the Appendix, we give some details of our numerical implementation of 
the frequency dependence. 


2. The fermionic many-body problem and the anharmonic oscillator 


We are interested in a 2PI fRG approach to fermionic quantum many body problems. 
The starting point for the definition of suitable generating functionals is the grand 
canonical partition function of a system of interacting fermions written in coherent 
state functional integral representation [47] . 




( 1 ) 


where the action is 

^['^] 2 ^ ^ '} Ota' ^ V 

otcx' (y.iQL'^a. 20 L 2 

Following Refs. |3S1 [35], we use a fully antisymmetrized notation here. This is 
customary in 2PI approaches and allows to derive the Bethe-Salpeter equation in 
arbitrary channels, compare Sect. The notation is based on indices a = {c,T,y) 
with charge index c = ±, imaginary time r and single-particle quantum number y. 
The fields '4>-,Ty = ''PyiT), 4’+,Ty = tS'ke Grassmann number values. The inverse 

free propagator is antisymmetric, and satisfies 

^(+,riyi)(-,T2y2) = ~ ^2)(9r2 + 


(3) 
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where e is the matrix of {Hq — fiN) in the single-particle basis. Here, Hq 
denotes the non-interacting part of the Hamiltonian H. The two-particle interaction 
vertex U is fully antisymmetric under permutation of its indices, Ua^a^a^ai = 

In the following, we study various 2PI fRG methods. Since the application of such 
schemes to many-body problems (like the Anderson impurity model) is involved, we 
here employ them to a simpler toy model in a first step. For that purpose, a “classical” 
choice Him 021113] is the quantum anharmonic oscillator with Hamiltonian 




-k 




9 .-4 

r ■ 


( 4 ) 


The position eigenstate path integral expression for its canonical partition function is 


Z = tre-^^= / D[a 




( 5 ) 


with real fields x, a suitably normalized measure D[x\, and the corresponding action 

(see e.g. Ref. [17]) 

l2 


S= dr- 


[drX{T)] 


--kH(x(r))| =J^ dr l^x(r) (Ordr+u^^ 2;(r)-k^ [a:(T)]'‘| . ( 6 ) 


Here, dr {dr) denotes a partial derivative that acts to the left (right). If we define the 
totally symmetric quantities 

( 7 ) 

( 8 ) 


C ^(r, t') = —drS{T — T')dr' — S{t — t')u 

C7(ri,T(,r2,T2) = 5{t[ - ti)<5(t 2 - ri)d(r2 - ri)g. 


'G> 


and identify r = a as well as x = tp, the action has formally the same structure as 
given in Eq. ([^ for the many-body problem. Due to the missing charge and state 
indices, computations for the quantum anharmonic oscillator are analytically simpler 
and numerically faster than for actual many-body problems. Typical observables can 
even be computed numerically exact which allows to assess the quality of approximate 
methods. Of course, the results cannot be directly transferred to many-body problems, 
since the fields are of different nature (real fields vs. Grassmann ones) and the 
additional degrees of freedom in many-body problems often lead to much richer 
physics. Nevertheless, basic properties of different approximation procedures can be 
tested on the anharmonic oscillator. 

The general derivations in the following sections hold for the fermionic many- 
body problem as well as for the anharmonic oscillator. Occasionally, the different 
commutation properties of Grassmann and real numbers lead to different sign 
prefactors. Then, ^ = —1 applies to fermions and (p = -1-1 to the anharmonic oscillator. 


3. Generating functionals in the 2PI formalism 

In this section, we briefly recall some basic relations for the generating functionals in 
the 2PI formalism and introduce the corresponding notation. We start by defining 
the functional 

Z[J] = = J . (9) 

The external source J is coupled to two fields here (cf. Refs. 07105]), which is decisive 
to obtain 2PI quantities by the Legendre transformation given below. J is defined to 
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have the same (anti-)symmetry as C~^, that is Ja'a = CJaa'- The first functional 
derivative of W[J] leads to the (time ordered) expectation value of ipa'^a’, that is to 
the propagator, 


= (V'aV’a')j = 

OJrvrv' 


( 10 ) 


In the following, it turns out to be convenient to use a combined index 7 = {a, a'). 
The higher derivatives 

( 11 ) 


= 


'■yi ■ ■ ■ 'J'Jyn 

are the 2n-point Green functions that are connected if each pair (cxj, a') is considered 
intrinsically connected. They obey the symmetry relations 


C[/(") ^ w(n) 

71---771 7P1---7P77 


w, 


(n) 


iai,a'J,72...7r 


= Cwt( 


Q:i,Q:ij.72---7n 


( 12 ) 


for arbitrary permutations P. 

The Legendre transformation of 1T[J] leads to a functional of the free variable G 
which is called 2PI effective action. 


r[G]= {-W[J]-J-G}\j^^^. (13) 

Here we used a dot-product notation, 

J ■ G = - J^Gj = - Jaa'Gaa' = ^ tl JG = G ■ J, (14) 

7 ctoc' 


where the factor ^ is related to the (anti-)symmetry of the factors under exchange of a 
and a'. As the quadratic part of the action is given by {G~^ -I- J) ■ {il’ip), it follows that 
r[G] = r[G, G] depends explicitly on G, namely due to an addend G“^ -G = | tr G~^G. 
In the non-interacting case, one finds 

Wo[J] = ln[det(-G-i - trln(-G-i - J), (15) 

ro[G] = [trln(-G) - tr(G-iG - 1)]. (16) 

The Luttinger-Ward functional is defined as 

$[G] =r[G]-ro[G] (17) 


and does not depend on G. The first functional derivative of T and is 

= ® = -t7G], (18) 

® = G-^ - C-^ - UG] = -E7G], (19) 

where E denotes the self-energy. In the last step, we used that the effective free (i.e. 
non-interacting) propagator in the model with external sources is given by (G“^-|-J)“^. 
As E[G] is the sum of all amputated one-particle 2PI diagrams with full propagator 
lines, it follows that 4)[G] is the sum of all vacuum 2PI diagrams (which are all 3PI) 
with full propagator lines. This expansion is shown in Fig. and given by 

7172 


1 

48 


E 

71727374 


Uaia2a3aiUa'^a'^a'^a'^Gj-^G^2G.y^Gj^ +0 (U^G^) . ( 20 ) 






+ 


+ 0{U^G^) 


Figure 1. Diagrammatic expansion of the Luttinger-Ward functional $[G]. 


This was derived e.g. in Ref. m for fermions and holds as well for the anharmonic 
oscillator. The individual addends are proportional to In this expansion, 

it is apparent as well that does not depend on C. The higher derivatives 
r("-) = S'^T/SG ... SG and = S^^/SG ... SG satisfy symmetry relations that are 
analogous to Eq. ( |l^ . 

Physical quantities are obtained at vanishing external source J = 0. We mark 
those quantities by overlining, hence J = 0,Z = Z[J],G = G[J], etc. For the Legendre 
transformed functions of the free variable G, the physical state is defined equivalently 
by = 0 and we write T = r[G],<I>*-^^ = and so on. For systems with 

spontaneously broken symmetry there may be several physical states. 


4. G-flow 


4-1. Flow parameter, flow equations, initial conditions 


In this section, we analyze 2PI fRG schemes where the flow parameter is introduced 
in the free propagator, G —>■ Ga (cf. Ref. [3S]). This method will be referred to as 
“G-flow”. As in IPI fRG schemes [T], we choose the A-dependence such that Gx, = 0 
at the beginning of the flow which ensures simple initial conditions for the flowing 
moments of $ (see below). We enforce that at the end of the RG flow Ga,. = G such 
that the cutoff-free system of interest is restored. 

The A-dependence of the free propagator Ga makes the action and the functionals 
depend on A: 5” —)■ Sx,Z —)■ Zx, etc. From now on, we suppress the index A most 
times to avoid overloading the notation. The flow of the functionals is given by 


Z[J]: 

W[J] 


sz - I 
SJ 

| = WW[J]-G- 



J[G] 


sw 

Jj 


. J - J. G = -W 


J[G] 


= GC- 


( 21 ) 

( 22 ) 

(23) 


Here, the dot denotes the derivative with respect to A, like Z = dZx/dX and 
C~^ = d{Cf^)/dX. We observe that P is independent of the interaction, P = Pg 
(recall that G is the free variable of P and thus independent of the interaction). 
Consequently, the Luttinger-Ward functional $ = P — Pg constitutes an invariant of 
the RG flow. 


$[G] = 0. 


(24) 


This is expected since 4> does not depend on G which carries the flow parameter [cf. 
Eq. (20)]. The flow equations for the physical values of the moments of $ are thus 


$1 


7l---7r 


= [(^1 = 1 
dX ^ 2 Z-^ 


^Gj, 


(25) 
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or briefly • G for n = 0,1, 2,..which is an inhnite coupled hierarchy 

comparable to those known from the IPI fRG [T]. 

The propagator G that enters the flow equations can be computed from G = 
(C-i-E)-i, 

G = -g(c-^ -t'jG, (26) 


a result that depends in turn on S = — 
equation for the self-energy, 


E = . Q 


(c--e) 


G. 


This leads to a self-consistent flow 


(27) 


An alternative approach to computing G is given by G = ■C~^ 

which is derived in Sect. 6.1 conmare Eq. ( [7^ . This, however, requires the solution 
of the Bethe-Salpeter equation (72) for 

We next discuss the initial conditions. From G = {C~^ — E)“^ and G\. = 0 


follows Gx- = 0. From Eq. (|20|), we deduce 

= 0 , 


<? = o, 


= U, 


= 0 , 


(28) 


71727374 


2 l-(i [^CKi<A:2a3a4 -f C(ai ^ ot'i)] 

+ C(«2 ^ a^)} + C(«3 ^ 013)) + C(q^4 ^ 04)] )(29) 


and more generally for n = 1 , 2 ,..., 

^C2n-l) ^ ^ ^C2n) ^ ^ 

Note that the nonvanishing result for ... contradicts Ref. [35] where it is 

claimed that = 0 for n > 3. 


4-2. Truncated flow equations and equivalence to self-consistent perturbation theory 

A straightforward truncation of the infinite hierarchy of flow equations for the 
leading to a closed set of equations is obtained by neglecting the flow of (and 

higher moments) for some fixed given m by setting As vanishes 

for odd TO, truncations at level to = 21 and to = 21 -|- 1 are equivalent, yielding both 
^ jji_ 2 PI fRG, this way of truncating is motivated from a weak 

coupling perspective: if U is sufficiently small, the impact of the higher moments of 
$ on the flow of the lower ones can be neglected. The relation of this truncation to 
established self-consistent perturbative approximation schemes is discussed next. 

The choices to = 0 and 1 are pointless since they lead to $ = 0, = —E = 0 

during all of the flow. Let us hence consider the case to = 2 or 3, which means setting 
= U. The resulting flow equation = —S = U ■ G can directly be integrated 
leading to T, = —U ■ G which is precisely the Hartree-Fock self-consistency equation 
for the self-energy. Integration of the flow equation $ = —E-G = G- ?7-Gin turn 
yields $ = ^G ■ U ■ G, which is the first order perturbation theory result for the 
Luttinger-Ward functional, compare Eq. ( |2^ . 

Diagrammatic approximations to the Luttinger-Ward functional combined with 
the self-consistency equation S = — are well known |33|. They lead to 
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conserving approximations [831 182]. In particular, the first order approximation to 
$ reproduces self-consistent Hartree-Fock, that is mean-field theory. More generally, 
the nth order approximation to $ leads to nth order self-consistent perturbation theory 
for the self-energy. This means that the self-energy is computed from diagrams up to 
order ?7”, with lines denoting full propagators G that are determined self-consistently. 

Now, we observe in general that truncating our set of flow equations at level 
m = 2n (or, equivalently, m = 2n-|-l) and solving the remaining flow equations results 
precisely in nth order self-consistent perturbation theory. This holds independently of 
which A-dependence is chosen for C, as long as C\. = 0 and Ca,. = C. In particular, 
this is independent of whether or not Cx regularizes diagrammatic expressions which 
are infrared divergent with bare propagators. Furthermore, this is true not only at Af 
but even during all of the flow (if perturbation theory is formulated with Cx instead 
of C). For the proof let us denote the nth order perturbation theory approximation to 
$[G] by A[G\. From A being independent of C follows A[G\ — 0 and • G. 

The initial conditions are 

, / = 0,1,..., 2n, = 0, Z > 2n -h I. (31) 

Obviously, the hierarchy of flow equations satisfied by the is finite and identical to 
the truncated exact hierarchy. As A and the truncated $ fulfill the same differential 
equations and initial conditions, they are the same. This concludes the proof. 

We thus derived a fundamental equivalence of two different approaches to the 
many-body problem: the solution of the 2PI G-flow equations truncated at level 2n 
gives exactly self-consistent perturbation theory to level n. It is rather remarkable that 
the resummation inherent to the RG procedure leads in the most natural truncation 
of the 2PI G-flow to a well known conserving approximation. This must be contrasted 
to most other truncated fRG schemes. E.g. in IPI fRG, the lowest order truncation 
for the self-energy is not equivalent to self-consistent Hartree-Fock [1]. Also, the other 
2PI fRG schemes discussed here (internally closed G-flow and two U -flow schemes, see 
below) are not equivalent to self-consistent perturbation theory. 


4..3. Internally closed set of flow equations 


Given the equivalence of truncated 2PI G-flow and self-consistent perturbation theory, 
we face the question whether other approximations to the G-flow are conceivable which 
lead beyond this known approximation to the quantum many-body problem. Such an 
approximation was studied in Ref. [35] . however not motivated by going beyond the 
above derived equivalence ||1] The idea is to truncate the hierarchy of flow equations for 
the by finding an approximate expression for in terms of The procedure 
is based on the leading (second order) perturbative contribution to 


(J)(3) 

717273 




2nd 


a2a30iA 


+ C(cn ^ 0^1)] 


74 


+ C(Q ^2 t-A 02 )} + C (®3 ^ 0 ^ 3 )) Gj^, (32) 

compare to Eq. (20). In Ref. [35], it is proposed to replace U in the expression 

for observe that this idea contains some ambiguity: while U satisfies the 


full index permutation (anti-)symmetry Ua^c 


= rc/< 


ap 1 Q;p2 CkPS CKP4 


, the symmetry 


II This equivalence is not observed in Ref. presumably due to the wrong initial conditions [see 
the sentence following our Eq. l |30[ l]. 
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of is only analogous to Eq. ( |l2| ). In Ref. [35], some specific order of indices of 
the interaction amplitudes U in Eq. (321 is chosen before replacing U —)■ Other 


choices would lead to other approximations. We consider it appropriate to replace 
instead 


U, 


Ol\0!.2C<C30la 


cj)(2),sym ^(2) , (j)(2) , (j)(2) 

Q:iQ: 2CK3CH4 ^ ^ 0:10:20:30:4 ' 0:20:30:10:4 ' 0:3040:20:4^ 


(33) 


where has the full index permutation (anti-)symmetry. The resulting flow 

equation, 

¥(2) = f$(3) ) _ (34) 

V 2nd/(7->.$(2).syn 

closes the hierarchy. Note that the (anti-)symmetrization of on the right-hand- 
side of the flow equation is indeed required since resulting from the flow equation 
does not have the full (anti-)symmetry. We refer to this scheme as “internally closed 
C-flow”. 

Due to = U the flow of ($^^^| 2 nd)( 7 ^¥< 
the flow equations then leads to 


starts at Integrating 


(j)( 2 ) _ (j)( 2 ) 

Af 


2nd 


Sa, = s 


2nd 




+ 0{U^G ), 
0{U^G^), 


(35) 

(36) 

(37) 


where 4>|2nd denotes the contributions to $ in first and second order in U, compare 


Eq. (20). Consequently, the solution of the internally closed C-flow comprises second 


order self-consistent perturbation theory. For the self-energy this means explicitly 

Sa, = Sexact + 0(173G2 SCPt)> (38) 

where G 2 SCPT denotes the full propagator of the physical state in second order self- 
consistent perturbation theory. 


GG Computing the thermodynamic potential 

One of the key observables depending on the interaction is the expectation value of the 
energy, which follows from the thermodynamic potential fl[G] = T[G]/j3 with physical 
value 




(39) 


For instance, the ground state energy is Egg = lim^_,.oo D. 

The thermodynamic potential can be computed from the values of $ and E 
that result at the end of the truncated or internally closed G-flow (or which follow 
in perturbation theory to $[G] from the self-consistency equation = _y; = 

). The physical value of the effective action is 


G ^-G-i 


r = $ -h ro[G] = 4> - ^ [tr In(-G) - tr{G-^G - 1)]. 


(40) 


With § tr(G-3G - 1) = (G-^ - G 




) • G = E • G, this means 


" = ^ 


-|trln(-G)-k$-f E. 


G 


(41) 
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For example in C-flow truncated at level m = 2, that is in lowest order 
perturbation theory to the resulting self-consistency equation E = —U-G 

corresponds to the mean-held approximation, G = Gmf- From $mf[Gmf] + Emf • 
Gmf = “^Gmf • U ■ Gmf follows 


Omf — 


— ^ tr ln(—G mf) — ^^mf ■ U ■ Gmf 
1 ^ 


= f^o[G]c._,GMF “ ^ 


(42) 

(43) 


which is the well known mean-held result. In order to evaluate one 

replaces the bare single-particle energies by mean-held dressed energies in the analytic 
result for no[G] = —{\nZo)/p. For the anharmonic oscillator, e.g., one replaces 

WG — > \J oJq + Emf in IuZq = — ln[2 sinh(/3a;G/2)]. 

Higher order truncations of the how produce a frequency dependent G which 
makes the evaluation of trln(—G) less trivial. Then, it is convenient to compute the 
thermodynamic potential from a how equation instead. Consider any approximation 
to the how equations for ... that result in some function 4)^^^ = —E^. This 

function allows to compute the how of G = (C~^ ~ ^ ^ 

G=-G(g-i- t) G, t=-S-G. (44) 


In the same approximation, the how of Q is 

Now fo[G] = G • G-i and rJ,^^[G] = G"! - G”^ = S result in 

n=^G-G-^ (46) 

We observe that the form of this how equation is independent of the approximation 
chosen for the how of E^. Hence, it is valid as well in case of the exact how, compare 

Ref. [55] . 

The initial condition 

Hy = 1 ($ + ro[G]),^ = -^[trln(-G) - tr(G-iG- I)],^ (47) 

is difficult to handle since Ga, = 0 leads to a divergence of ln(—G). Following Ref. [35], 
we study instead 

An = n-Ho[G] =n-k^trln(-G) (48) 

Ip 


with initial condition AHaj = 0 and how equation 

An = ^ An = n + tri„(-c) = i(G - c) • c-. (49) 

At the end of the how, we recover the thermodynamic potential as Haj = Hg [G] | -|- 

A^Af • 
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^.5. Application to the anharmonic oscillator 


We test the performance of C-flow truncated at level m = 2,4 or internally closed on 
the quantum anharmonic oscillator (with C = + a = r). Due to translational 
invariance in imaginary time, it is preferable to work in frequency representation. 
Hence, we Fourier transform the m-point functions to bosonic Matsubara frequencies 




— ILOn-. Ti ...—Tr 


$ 


nin-^ 

(m) 


(ti, .. .,T^) 




where dr. Mind the choice of signs in the exponents. Here, 

is a representative of propagators (including G), whereas represents vertices 

(including E). For the free inverse propagator one obtains 


)nn' — /ddn+n',0 {i^nl^n' 4” ^g) — /dSji+n',0 (^n 4~ ^g) ' 

More generally, all relevant functions are proportional to j36ni+...,Q- We refer to the 
remaining part by skipping the last index, e.g. 

^nn’ = Pdn+n',0^n, ~ 


For the flow parameter dependence of C, we choose specifically 

= ~ {‘^n + Wq + a) , 


(54) 


with Ai = 00 and Af = 0. This can be considered as an analogue to the so-called 
hybridization flow parameter for single-level quantum impurity problems [laiiiiiii], 
which has the form = [ja;„ — e-I-* sgn(a;„)(r-|-A)]~ for a system in grand 
canonical equilibrium. Here, e denotes the on-site energy of the level and F the 
hybridization to the leads. Generally speaking, the results following from the solution 
of truncated fRG flow equations depend on the cutoff choice (in contrast to the solution 
of the infinite hierarchy of equations) [T]. Therefore, it is an interesting observation 
that the results of truncated G-flow are completely independent of the particular 
choice of the flow parameter as they are equivalent to self-consistent perturbation 
theory. This does not hold for the internally closed C-flow. By comparison to exact 
results, it was shown for selected fermionic many-body models that the hybridization 
flow parameter is a good choice mus] in IPI fRG. This provides the motivation 
to use its analogue for the problem at hand. We here do not further investigate the 
influence of the cutoff choice on the quality of the results. 

Dyson’s equation yields 


Gn — [(C ^)n — — — (iOn + ^G + -^ + 


(55) 


Due to C„ ^ = — 1 the flow equation (49) for the thermodynamic potential takes the 
form 


1 


^^ = -^T.iGn-Cu). 

^ n 

The flow of the self-energy E = —reads 

f __C 


(56) 


(57) 
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with 


= -G„ 




Gr, 


(58) 


The initial conditions are = 0 and < 1 )*'^^’^' = q. In case of G-flow truncated at 

^ nin-^n2 ^ 

level TO = 4, the solution is self-consistent second order perturbation theory, 

5 ,5^ 


^"1 - 2 / 3 ^ "" fi/12 


6/32 


^ ^ Gn'zGn^G 

n2’T.3 


(59) 


This equation is solved self-consistently for each A (which is numerically faster and 
more accurate than integrating the flow of S, and $ 1 ^ 1 ) and the result is used to 
calculate the flow of All. 


$ 


For the internally closed G-flow, Eq. (571 is solved numerically together with 

=_lvfe) 


( 2 ) 

‘ n\n\n2 


n3 


ssym ^ 

"ni,n2,n3^ 


+ (*21 ^ ^^l)} 


1722 =—ni—rij^—n2 


G„3,(60) 


^(2).sym „ _ 

^ni,n2,n3 g 


-p $( 2 ) 
"l.''l2,"3 »l2,"3,’2l 


. $( 2 ) 
^n3,ni,n2 


Eq. (601 exploits the symmetry = E _„3 which implies G„j = G-r 


(61) 
It follows 


—( 2 ) —( 2 ) 

from Eq. (12), which translates to n 2 ~ n'nin 2 in frequency space, via Eq. (57). 

We restrict our numerical investigations to the case of vanishing temperature and 
thus take the limit /3 —>■ 00 . The Matsubara frequencies become continuous then, 
uin —>■ and sums turn into integrals, / If- Transforming, for instance, 

Eq. (59), one finds explicitly 

) 7 

G{v2)G{vz)G{-v-V2-V3,)- (62) 

.... ..... „o 27r 

We here focus on two observables. One is the shift in the ground state energy 
induced by the anharmonicity, Cgs = Egg—E^^ = lim^_>oo Ail, with the non-interacting 
{g = 0) ground state energy Eg^. The other is the fluctuation of the position g = 
which corresponds to the particle density in a fermionic many-body problem. In the 
following we therefore refer to g as the density. At the end of the flow, it can be 
extracted from the self-energy as 






' diy2 
^ 2 tt 


0 = -G(O,O) = 


f 


dv 


“len dl/ 


1 


1 


27r i 


-Wr 




TT -|- Wq -|- 


(63) 


Here, the last integral refers to the numerical evaluation on a grid of frequencies v„ 


with n = —TOien, ■ • ■, "mien and r'„ = —v-n (details on the grid are given in Appendix 
0 . The integral is then calculated by an interpolator. The very last addend is an 
approximation to the missing integral part — ^G{v). 

In case of G-flow truncated at level to = 2, the solution is mean-field theory. At 
zero temperature, the self-consistency equation Emf = —U ■ Gmf takes the form 


V 5'/ 2\ 9 

^MF “2'* ImF — ^ 


I 


(64) 


.^MF 


In this case, we do not integrate numerically any flow but solve the self-consistency 
equation directly at Af and compute the thermodynamic potential from Eq. (43) 
00 ). This leads to the regular mean-field result [4] 

I 


(taking (3 


e^^ = - 


COr 


Smf — 

25 


WG 


MF ' 


(65) 
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—( 2 ) 

For the internally closed C-flow, ’ depends on three frequencies. Consequently, 
the number of flow equations scales ~ ™fen- sum must be performed on the 


right-hand-side of Eq. (60), the effort of computing all the right-hand-sides of the set 
of differential equations in each step of the flow scales ^ This severely limits the 
accessible mien in this scheme. In comparison, the regular C-flow schemes truncated 
at arbitrary order m always produce ^ mien number of flow equations. This is because 
only An and flow in the regular schemes as the flow of can be determined 
self-consistently. Nevertheless, there is (of course) an increasing cost for calculating 
increasing orders m: the effort to calculate the right-hand-side of the self-consistency 
equation for Y,{v) increases. For example at m = 4, the effort to calculate all equations 
for one step of the flow is ~ mfg„—this can be managed by pre-calculating one of the 


dv 21 (ir' 3 -integrals in Eq. (62) before actually calculating the right-hand-side of the 
equation itself. 


4-6. Numerical results for the anharmonic oscillator 

For all numerical considerations, ujq is set to 1 in this paper. Fig. shows numerical 
results for the ground state energy shift and the density as functions of the coupling 
constant g. The data obtained from C-flow truncated at levels m = 2,4 or internally 
closed are compared to the exact results. As the actual (“absolute”) curves are very 
close to each other, the main plots show the relative differences to the exact result. 

We compute the numerically exact data by matrix diagonalization: the full 
Hamiltonian (with quartic anharmonicity) can be expressed as a matrix in the basis 
of the energy eigenstates |E°)) of the unperturbed harmonic oscillator. Diagonalizing 
this matrix yields the energy eigenvalues and eigenstates of the anharmonic system. 
The ground state density is then obtained as = {Eo \ |Eo), where the matrix 

(^ )m,n 

elements: 


= {E°m\ I A'n) in til® unperturbed eigenbasis only has the following non-zero 
(i^)n,n = n + ^ and (x^)„+2.n = (i^)n,n+2 = ^y/n^y/rT- 


= n + ^ and (x^)„+ 2 ,n = {x'^)n,n +2 = + ly/n + 2. This can 

easily be proven from the ladder operator representation x = (^a + /v^- A matrix 

size of 200 x 200 is found to be sufficient for converged numerics. 

In all numerical C-flow computations, the frequency grid covers the range from 


= —5120 to t'mien = 5120 with 2mien +1 grid points, see Appendix A for details 


on the grid. The truncated C-flow is computed with mien = 500, while the numerically 
much more demanding internally closed C-flow is solved with mien = 45. In the latter 
case, comparison to mien = 35 and 40 data suggests that convergence was reached 
on the scale of the figure^] This is also supported by the observation that data from 
truncated C-flow computed at mien = 45 are close to their converged counter-parts. 
For all C-flow methods, the flow starts at A"™ = 10^®. Convergence with respect to 
Afum requires such a large value since the ground state energy depends on Af“™ in 
leading order via an addend proportional to l/\/Af“™. 

Obviously, the self-consistent Hartree-Fock approximation, that is C-flow 
truncated at m = 2, yields already quite good results with a relative discrepancy 
of less than 3% for g up to 50. For the problem at hand, this is plausible as mean-field 
theory can be expected to capture reasonably well the effect of the cc^-anharmonicity 
on low-energy properties by self-consistently optimizing the frequency of the harmonic 
oscillator. The improvement achieved by second order self-consistent perturbation 
theory, that is C-flow truncated at m = 4, is rather low in comparison: the data is 


% If at all, the energy curve is expected to move a little closer towards the exact data for higher mien- 
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Figure 2. Numerical results for the ground state energy Cgs (a) and the density 
Q (b) as a function of the strength of the anharmonicity g. The lines show data 
obtained from C-flow truncated at level m = 2 (dashed), that is mean-field theory 
(MF), from C-fiow truncated at level m = 4 (dashed-dotted), that is second order 
self-consistent perturbation theory (2SCPT), and from internally closed C-flow 
(dotted; ICCF). The main plots show the relative difference to the exact result, 
whereas the insets show the absolute mean-field and exact curves. 


typically 1 to 2% off for Cgs, while g is reproduced even less accurate than by the 
Hartree-Fock approximation for g > 20. The results from internally closed C-flow are 
significantly better, with relative deviations below 1% for both observables. The idea 
of internally closing the set of differential equations by expressing in terms of 
is thus the most successful C-flow scheme investigated here. However, the price to pay 
is a large computational effort ^ ^fen computation of the set of differential 

equations in each step of the flow. Additionally, it is unclear whether this approach 
can be systematically generalized to higher orders. Therefore, it is not possible to 
compare results from truncations at different orders in order to assess the quality of a 
given approximation, a procedure often used in fRG studies [T]. Anyhow, hypothetical 
higher orders of internally closed C-flow seem out of reach for numerical computations 
as the scheme implemented here is already numerically quite demanding. 

The question emerges whether one can think of an RG scheme that better fits 
the 2PI nature of the general approach. One possible answer is the C-flow that shall 
be discussed in Sect.|^ Sect, [^provides some preliminaries for this discussion. 


5. Pair propagator and Bethe-Salpeter equation 


In the discussion of the C-flow scheme below, some notions related to four-point 
functions and the dot-product introduced above will be useful. They are provided in 
this section which can be considered as a sequel to Sect.[^ 

From the non-interacting functional Wo[>-^] = —2 ti'ln(—C“^ — J) follows Go[J] = 
= (C-i -k J)-^ From = -SGo[J]/SJ follows [Jo[C]] = n[C], 

with the pair propagator 


<^^72 _ ^ 
- 1 


^G-\, 

OL\OL-^ 


— c , c 

^cx.-^OL2 ^a.icx.2 


CGo:ia2 Ga' aL 


^ 71 72 [^] 


( 66 ) 
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In case of a particle number conserving system, 11 denotes either a particle-hole or a 
particle-particle propagator, depending on the charge indices. For example the charge 
indices ci = —c'^ = C 2 = —c '2 render a particle-hole propagator, while the charge 

indices ci = = —C 2 = —c '2 lead to a particle-particle propagator. 

In the general interacting case, one has 


where 


= 


r _ ^'A'2 


6J 


5G SJ 

SJ ^ 

(67) 

SG SJ 6G 

II 

1^ 

I'd 

C^aia2^a'^a2 


(68) 


is the neutral element with respect to the dot-product, X ■ I = I■ X = X. We conclude 


p^(2) ^ (69) 

and in particular Fg ^ = 11“''. Here, we introduced the dot-product inverse of four- 
point functions that satisfies 


An example is 


(X ■ i V X 

^ ''7172 9 


(n“'') [Gl = g7 g 

V / -71 -yo L J CKi CXo < 


^inv _ r 

7173^7372 ~ -*7172 ■ 

(70) 

aia '2 ~ ^ 

(71) 


where G ^ denotes the regular inverse with respect to the a-indices. Now, F^^^ = 
$( 2 ) -(- Fg^^ = -I- H'"'' leads to the Bethe-Salpeter equation 


p{/(2) ^ (Jjinv ^(2)ynv ^ . ^(2) . n ± . . . = n - n • . P17(2) 


(72) 


6. G-flow 

6.1. Flow parameter, flow equations, and initial conditions for the plain U-flow 

In IPI fRG, the flow parameter is typically introduced in the one-particle propagator. 
In G-flow, this procedure was carried over to the 2PI case. We next explore the 
consequences of introducing the flow parameter into the interaction, U —> Ux, as 
done in Ref. [3S]. This method will be denoted as “G-flow”. We set Ga, = 0 at the 
beginning of the flow, such that the flow starts at the non-interacting model. At the 
end of the flow, U\^ = G restores the original interacting problem. We maintain the 
full index permutation (anti-)symmetry of the interaction during all of the flow, such 

that bGiQ2Ct3Ct4 — C 

In Ref. [7], an “interaction flow” method in the framework of IPI fRG was studied. 
The authors describe that a flow parameter G —>■ AG can be substituted to a flow 
parameter G —>■ A^G by rescaling the fields. One could jump to the conclusion 
that there is in principle no difference between G-flow and G-flow, at least for a 
multiplicative flow parameter. The rescaling of the helds, however, leads to a A- 
dependent rescaling of all Green and vertex functions. Thus, the structure of the flow 
equations for the non-rescaled vertex functions differs for the two flow parameters and 
suggests different truncation procedures. 

It is conceivable to combine the G-flow and G-flow schemes by introducing a flow 
parameter into both, G and G, compare Ref. [35] and also [15]. Here, we discuss only 
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the pure t/-flow. In a combined scheme, the flow equations will include additionally 
the terms derived for the C-flow in Sect. 14.ll 

Possible infrared divergencies of perturbation theory ought to be regularized by 
the chosen flow parameter. This can indeed be achieved in the framework of the U- 
flow. It has been shown e.g. in the context of partial bosonization for one-dimensional 
models with interactions that involve only small momentum transfers that a cutoff in 
this transfer regularizes all infrared divergencies [50] (see also [36]b 

Ux induces a A-dependence of the action and of the functionals. We now derive 
the flow equations for the . We refer to this scheme as “plain C/-flow” in order 
to distinguish it from the modified scheme that is proposed in Ref. |36j and which we 
describe in Sect. |6.3| The flow of the functionals is given by 




1 

3! 


Trt/ • 


S^Z 

5J5J' 




. t/. tT(i) 


f 


= -^l.7[G] = 


3! L 


Tr U 




+ G-U-G 


(73) 

(74) 

(75) 


In the last expression, IT^^^ [J[G]] = (r(^^[G])’”''. We used the “double-index” trace 


1 


(76) 


compare Ref. |30|, which is to be distinguished from the “single-index” trace trT = 
Yaa- As the non-interacting effective action Tq does not depend on U, the flow of 
the Luttinger-Ward functional $ = T — Tq is given by 

$ = f = iTrt/- + , (77) 


where W^^'> = -h 

In order to compute the flow of the given by = $(") need 

expressions for G and For G, we use Eq. (26) for the practical application to the 
anharmonic oscillator. However, for gaining general analytical insights, a different but 
equivalent expression turns out to be useful: from 0 = F*^^) = F^^^ -|- F’^^) . q follows 
G = . f(i) = . $(i). (78) 

For = S^/SG, we find 


$(i) = ■ 

3! 




/(5W(2) 


V 


6G 

tF( 2 ) . 


I (5n 
26G 


H” 


SU 

SG 


n™'' - $(3)). 14/(2) 


I SU 

2Jg 


(79) 

(80) 


writing indices explicitly, this means 

4>W = ^ XI 1 n”'' • U ■ H'"'' 

^ 72 V 



G^ 


CtiCt2Gt2^1 


(81) 
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(82) 


This yields explicit expressions for G = — and 

$ = <!)-(- $(1) . ([j = (j) _ 

as well as 

^(1) = $(i) +$(2) . 4,(1) ^ (^/ + $(2) . ny'"' • (83) 


In order to evaluate the right-hand-side of Eqs. (82) and (83), we need to know 


[which also allows to compute = (ll*”'" -|- and [cf. Eq. (81)]. 

This is the onset of an infinite hierarchy of coupled flow equations, where $ depends 
on ..., and where depends on ..., <I)("+2) for n > 1. 

The starting point of the flow is given by Ux^ = 0 and corresponds to the non¬ 
interacting model. Consequently, = 0 and 


= 0 , 


c&W =0, n = l,2,3, 


The flow of the thermodynamic potential is given by 

n = t = 4 (r -f S ^ 


P /3 


/3 P 


(84) 


(85) 


It starts at the non-interacting potential, Vlx. = Oo[C'] = —4 InZp. 


6.2. Truncation of the plain U-flow 

In the plain t7-flow, the lowest truncation of the hierarchy of flow equations that leads 
to a flowing self-energy is at level m = 3: one sets 4)^"^ = = 0, that is to its 

initial value, for n > 3. The resulting flow equation -|- • G = reads 


$(2) = -Trt7 


2 TT(^) . n'"'' 

-h . n” 




SG 

<5G^ 


^ • n”" • 
SG 


n“''.w( 2 ) + 


1 

2 


( 86 ) 


and produces a with explicit time or frequency dependence that enters the flow 
of the self-energy. Solving the flow equations is then already quite involved. In order 
to find a simpler approximation, we use the expansion =11-1-0 and get 

s^Ti 


$(2) = _ Tr G ■ 
4 

1 ( 2 ) 


SG^ 


O 




= u 


l-b0 




(87) 


This can be integrated to -|- 0{U‘l) which we approximate by = Ux. 

This result can be understood also as the leading perturbative approximation to $(2), 


compare Eq. (20). The flow equations for the self-energy and the thermodynamic 


potential are then 


= -« E(^ + f^>-nA) 


inv 
7i 72 


72,73 


(I + Ux- Tlxr^ -Ux-il + Tlx- + 


Ux 

2 




( 88 ) 


a20C3Cx.^a2 
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(nr + C/a) 


C/aI +!^ 


(89) 


Let us address the question how the solutions of these approximate flow equations 
compare to perturbation theory. It turns out that they comprise first and second order 
perturbation theory (with bare propagators), 

Sa, = Sexact + O (C/3) , Ha, = Hexact + O (C/^) . (90) 


For the proof, we start from the initial condition Sai = 0, G\ = C, 11 a, = Ho. From 
t = 0{U)[l + 0{U)] follows E = C>(C/), hence G = C+ C>(C/),n = Ho+ C>(C/). This 
leads to 

= -c -U ■ (CEC) + C/ • Ho • t/ • G]^^ 

+ ^ • Ho • G + C7 • Ho • G)„,„3„A„iG^3 + 0{U^U). (91) 

73 


The first addend can be integrated to {—U ■ C) which is first order perturbation theory 
for E. We insert E = —U ■ C + 0{U'^) into the second addend, 

-U ■{CY.C) = U-Tio-U-C+ 0{U^U). (92) 


Together with the third addend, this can be integrated to the value of the mean- 
field-like (non-2PI) second order self-energy diagram C/ • Hq • G • G, plus third order 
corrections. The fourth addend can be integrated to 

g 'y (93) 

73 

which is the value of the 2PI second order self-energy diagram. Thus, we recovered 
all diagrams up to second order, and the remaining 0{U'^U) leads to some correction 
of 0{U^). The proof for the thermodynamic potential is similar. 


6.3. Modified U-flow starting at mean-field theory 


In Ref. |36j, a slight modification of the G-flow is proposed which results in a different 
starting point of the flow, namely mean-field theory instead of the non-interacting 
problem. The modified G-flow is defined by 


$A(modified G-flow) = $A(plain G-flow) -|- • (G — U\) ■ G, (94) 


where <I>A(plain G-flow) refers to the scheme described in Sect. 6.1 For an 


interpretation of the modified G-flow, consider the expansion of the Luttinger-Ward 
functional given in Eq. (20). In the plain G-flow, each G in this expansion is replaced 
by Ga. In the modified G-flow, the first order contribution 


• G • G is excluded 

from this replacement. At the starting point Ga, = 0 of the flow, ^a, in the modified 
scheme is given by precisely that first order contribution, 


$A = ^G-G-G, 


= -Ea. =G-G, 


= G, 


= 0, n > 3.(95) 


As E = — G-G is the Hartree-Fock self-consistency equation, we see that the flow starts 
at mean-field theory (cf. [27]). At the end point Ga, = G of the flow, the full genuine 
Luttinger-Ward functional is recovered. If required, the flow of the other functionals 
is derived from the flow of <I>a. The effective action is given by Fa = d)A -f Fq, and Wa 
is obtained from Fa by a Legendre transformation. 
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In the same spirit, one could define more general flow schemes that start at nth 
order self-consistent perturbation theory. In the definition of the replacement 
U —>■ U\ would be applied only to contributions of order 17"+^ and higher. Depending 
on the problem of interest, starting the flow with nth order self-consistent perturbation 
theory might or might not be advantageous; see the discussion of truncated C-flow for 
the anharmonic oscillator above and Ref. [51] for the single-impurity Anderson model. 

In the framework of the modified scheme, the functional <I>(plain [/-flow) does not 
play the role of the Luttinger-Ward functional. Therefore, we label it from now on 
such that Eq. (941 takes the form <1);^ = ■ {U — U\) ■ G. For the flow equation 

U ■ G = — \ Tt: U -II. We insert 4' from Eq. (77) and obtain 


follows $ = 4' — /G 


1 


4> = - Tr [/ • 
3! 


(rf'"' -k 4'(2)y 


-n 


(96) 


(l + 4'(2) . ny"""- [7 • (/ -f n • u\ G^^ + C>(4'(3)) , (97) 

^2 ^ Q:iQ;2CK2^1 


where 4^^^^ = 4?^^^ + U\ — U and 4'^^^ = On the analogy of Eqs. (82) and (83) 


follows 

i = $ - $(i). (98) 

As for the plain [/-flow, the flow equation for the thermodynamic potential has 
the form O = 4>//3. The initial value is now the mean-field potential Omf indicated 


in Eq. (43). 


6.4- Truncation of the modified U-flow 


In the modified [/-flow, the lowest truncation of the hierarchy of flow equations that 
leads to a flowing self-energy Ea = —"hA^^ is at level m = 2: set = U and 

= 4>^"^ = 0, n > 3. The flow of the self-energy in this approximation is 

^7. = -^E(^ + t^-nA)‘T7. 

72,73 


X 


(/ + [/a • Ha)'"" • [/a • (7 + Da • Ga)”" 


Ux 


Ot20c30c^Ct2 


^73’ 


(99) 


where U\ refers to the flow parameter dressed value of the interaction, while U denotes 
the bare one. The flow equation for the thermodynamic potential = 4>//3 reads 


f^A = 3^TrGA 


/—inv \ — 

(Ra + Ux) - Ra 


( 100 ) 


It can be shown that the solutions of these approximate flow equations comprise 
all diagrams contributing to self-consistent Hartree-Fock (mean-field) and to second 
order perturbation theory with mean-field propagators, 

^Af = Eexact+G , IlAf = Rexact + G ^G^G]y[p) .(101) 


For the proof, we start from the initial condition Ea, = Emf, Ga, = Gmf, Rai = Rmf- 
From t = 0{UUG^) follows Ea = Emf+G(G^gE), hence Ga = Gmf+G(G^gE), 
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Ha = Hmf + 0(C/^Gmf)- Therefore, 

^7i “ o 'Ux + Ux' Hmf • Ga) Gmf.73 + O (u^UGy^-p) ■ (102) 

O V / OL\OL20L'r,OL\ ^ / 

73 

Integration from Ai to Af yields 

£a, = Smf+^ E ■ nMF • Gmf,73+C1 (g'G^) -(103) 

73 

The second addend on the right-hand-side is indeed the 2PI second order self-energy 
diagram with mean-field propagators. The proof for fl can be performed analogously. 


6.5. Application to the anharmonic oscillator 


For a practical application of plain and modified t/-flow, we use Eq. (26) for the 
calculation of G. In plain [/-flow, one obtains for the self-energy (still for general 
7 -indices) 


7 = - E 


72 


inv Ux 


(I + Ux-Tlx) -Ux-il + Tlx-Ux) + 


Gt 


a.lCii2<y.20ii 


9 Ha.7172^^2 . 


(104) 


Introducing the quantity Ta = I — (/ -\-Ux • IIa )*'''', one finds 


7=- E 


72 


T A • [>A • TI - T A • [/a - Ha • TI -f 5 [/a 


H 72 + 2 

a\cx.20i'20L’-^ 






^ (105) 

(106) 


On the other hand, one finds in the modified [/-flow that 


H73 —3 


^ [Ta-[/a-TI-Ta-[/a-[/a-TI 


72 


, ,<-bE^7372G4, (107) 

Q:iQ; 2 Ck 2^1 ^ 


72 




(108) 


Note that due to (T]^).y 2 _.y 2 = Ta; 72 , 7 i, all occurrences of can generally be replaced 
by T. This shall be done below. 

Now, we again consider the specific problem of the anharmonic oscillator [i.e. 
C = -1-1, 7 = (r, r')]. The above equations shall be Fourier transformed according to 

Eqs. (50) and (51). Note that is neither propagator nor vertex but rather ~ I. 
Suitable Fourier transformation definitions are 

T 


- n\n\^n2n2 




/ ... / e 1 "1 ' = "2 =T(ri,ri,r 2 ,T 2 ), 

J n r 2 

... e 1 ^ "2 2T (Ti,Ti,r 2 ,T 2 ). 


(109) 


( 110 ) 


Let us specify the flow parameter dependence of Ux via a function /a S K with 
/aj = 0 and /aj = 1 such that 

HA(Tl,r(,T2,T2) = fxgSin =t{=T 2 = T^). 


(Ill) 
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In practical calculations, f\ = X with A, = 0 and Af = 1 is used. The choice of U\ 
makes it possible to show that 

'^X;nin[n2= 1 + X! -"'l -n'l > (1^^) 

^ n 

r 1 

%n,n',n2=h9GtGl, 1 + ^ . (113) 

As expected from one finds '^\-n 2 ,n,+n{-n 2 ,n^ = \-,n^,n'^,n 2 - 

From now on, we will only use T>,. It turns out that the third index of is always 
summed over independently. This sum depends only on the sum of the first to indices. 
Thus, we introduce the definitions 

= + '/a9(GG)„_„,, (114) 

(^)iV =-oY.G-nGn-N. (115) 

^ n 

With these findings, the flow equations in plain {7-flow become 


Oa = —- 


E G^,-„,G^^Ta;. 


Ag I A 
8 1/3^ ” 


Aa _ /Ag .f\9sr^^ T^A /Ag \^^A /AgM^^bA ^ 1 

12 o / A A;Ti^K-rii+ 0^ / . ^ X\n^n-m 2jj / / " 2jj / / 


The last addend in the second equation depends on and makes the equation a 
self-consistent one. However, this self-consistency problem can be solved explicitly as 
the addend does not depend on ni (i.e. it is constant with respect to the external 
frequency). With the ansatz + G^ where denotes the first three 

addends and the last addend, one finds: 

G" = - 1 + #E(GA)' ( 118 ) 

^2 J n2 

With the same ansatz E^ = -|- C^, one finds for the modified {/-flow: 

E ( 119 ) 

711712 

TVA _ fx9 'Y'2 , hgsrxi' T^A non\ 

^ni 12/? ^ ^ ^ A;7i^n—Til ' no / ^ ^ X;n^ri—ni'! 

^ n ^ n 

- fi + §E §E (121) 

^ n2 J ^ n2 

In both schemes, a non-trivial frequency dependence is generated in the self-energy in 
these lowest-order truncations (typically, this is different in IPI fRG schemes). Note 
that the following symmetries hold: G„ = G_„, (GG)^ = (GG)_^, Tn = T-n and 
Dn = G_„. 
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Figure 3. Numerical results for (a) the ground state energy and (b) the density 
as a function of the anharmonicity g. The relative difference to the exact results 
is shown. The lines present results obtained with the plain f/-fiow (dotted; PUF) 
and the modified U-fiow (light [red] dashed; MUF). The C/-flow data is calculated 
with = 180. The mean-field (MF) data are given for comparison. 


Now, the P —>■ oo (i.e. T —)■ 0) limit can be performed. Then, the initial conditions 
in plain 17-flow are 

= ^, rA.M = o. (122) 

The initial conditions in modified U-flow are 

^Ai = t/^MF ~ = Emf = t— , ■ (123) 

^9 4 + 

6.6. Numerical results for the anharmonic oscillator 

We implemented the zero temperature U-flow equations making use of the non-uniform 
frequency grid introduced in [Appendix A[ In both schemes, the number of equations 
scales ^ mien and the effort for computing the right-hand-side is also ~ mien; the 
effort for the calculation of the set of differential equations in each step of the flow is 
thus ~ n^ien (like in regular C-flow at truncation order m = 4). Reaching numerical 
convergence with respect to mien is thus not an issue as it can be tuned to rather large 
values. Fig. shows converged numerical data for mien = 180. 

The modified U -flow curve is a lot closer to the exact result than the plain U -flow 
curve. In fact, it reproduces the exact results within 1% and can thus compete with the 
internally closed C-flow method. Recall, however, that the effort to calculate the set 
of differential equations in each step of the flow scales ~ mf^j^ in the modified U -flow 
while it scales ~ fnfen i^ i-h® internally closed C-flow. There is an intuitive reason for 
the success of the modified l/-flow: the initial condition corresponds to the Hartree- 
Fock approximation which reproduces the exact results quite well already (< 3%). 
Thus, much of the interacting physics is captured initially and the flow only has to 
take care of the missing part. Apparently, starting the flow from the non-interacting 
system is a significantly more unfavourable initial condition which explains the less 
successful results of plain U-flow. This difference between the {/-flow schemes can also 
be understood from the respective discussions of the degree of perturbation theory 
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contained: Both schemes comprise second order perturbation theory but plain U-flow 
does so with the bare propagator while modified tJ-flow does so with the mean-field 
propagator. 


7. DFT-fRG 


A method similar to 17-fiow was proposed in Refs. [SI |S] and applied to the 
quantum anharmonic oscillator in Ref. |46| . Instead of working with an external source 
depending on two indices [like in Eq. ([^], a purely diagonal one = Saia^Jai 

is introduced. Thus, the formal starting point of the considerations is not a 2PI 
functional but rather a two-particle point-irreducible (“2PPr’) one. This term is 
introduced in Ref. [52] with the following descriptive explanation: a 2PPI diagram is 
a IPI diagram that cannot be split into two by cutting two internal lines attached to 
the same vertex. This implies differences in the derivation of the flow equations and in 
their structure. In particular, it turns out that the Luttinger-Ward functional cannot 
be computed reasonably because the relation G = — ^ cannot be solved explicitly 
for J any more. The manifest remaining option is to work with the functional P and 
its moments. Furthermore, the diagonality of J causes G to be the density of the 
physical system. According to the discussion of Sect. |4.4| the ground state energy 
can be expressed through the functional r[G]. Therefore, one works with an energy 
functional that depends on the density. This brings about a close relation to density 
functional theory [53] (in the plain Hohenberg-Kohn sense) and motivates us to use the 
term DFT-fRG. The physical state indicated by a bar is defined as above, namely 
by the minimum of the energy functional. Thus, this method yields the minimum 
of the energy and the corresponding density of the interacting system at the end of 
the flow. In this sense, DFT-fRG is a way of conducting density functional theory 
calculations alternative to the Kohn-Sham idea [53]. Truncated DFT-fRG can even 
be used to derive an approximation to the energy functional r[G] of the interacting 
system. It remains to be seen if this approximation—which has the structure of a 
Taylor expansion around the ground state density, see below—leads to useful results 
for fermionic many-body problems. 

Introducing the flow parameter simply via U —>■ XU, Ai = 0 and Af = 1, flow 
equations can be derived for P^") without further conceptual complications. However, 
the initial conditions T^"^ are rather lengthy to compute because even for the non¬ 
interacting case r takes a non-trivial form (in contrast to $). The explicit calculations 
can be found in Ref. [IS]; note that in that work the moments T^") are not defined as 
functional derivatives but rather only implicitly [namely via Eq. (14) of that work]. 
Thus, the P^"^ cannot be expected to obey all the index permutation symmetries that 
functional derivatives would haverd 


Once the flow equations P^"^ 


and the initial conditions Pl"'^ are determined, a 


natural truncation is obtained by setting T^”^ = P^"^ for all n > m for a certain m. 
This is again motivated from a weak coupling perspective (cf. Sect. 4.21. At the end 
of the flow, one obtains an approximation to the interacting P^"^ for n < m. Along 
with the non-interacting P)^ for n > m, these moments constitute an approximation 
to the energy functional which is valid “locally” around its minimum. In this way, 
DFT-fRG is able to provide an approximation to the interacting energy functional. 


■*" This implicit definition was dealt with consistently in Ref. SS] but the reader should keep it in 
mind in order to avoid confusion. 
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We here specify the (first three) flow equations for the quantum anharmonic 
oscillator. For a self-contained presentation we stick to the notation introduced above 
up to a redefinition of G\ including a factor ^ (compare to Ref. [15]). From F^^^ = 0, 
one obtains a flow equation for G\. Instead of providing a flow equation for F)^ , 

It Is equivalently given for W\ . Adopting the definition for F by functional 
derivatives and additionally performing the zero temperature limit (this implies some 
minor differences to the equations given in Ref. [15]), one finds 


= — 
^ 24 


poo 1 

dl'lTT7(2), 




Gx = (0) 


-Gx+ ^ 


27r 






(124) 

(125) 


V) = vr^ 1 + -GAr^^^(0, ix) - I FH^i-^1, ^) 
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12 

* '^^1 tI7(2V 


12 - 




^tI7(2V„. v(4V 
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wf(zzi+^) (rf(j^i,i^))'+iRf(i^i-i^) (rf(^^i,i^))^ 


The initial conditions are 
WG 


flo — 


2 ’ 


1 


= — 


1 


Tg \lXl, 1X2) — “Wq {yi + + V2 + 12 Wg) j 


• (126) 

(127) 

(128) 


rA'‘\l^l,J^2T3) = 


2wn 


V\-\-V2') -f (iZi-|-^'2) ^'3-|-^'3-|-12 wq j +1^11^2+ 1^2 “b 12 wq) 


[vi -|- V2 ) + 4ujq 




+ ^2) + 4a;Q^ ^(ia 2 + 1A3 ) -|- 4wq^ 


+ {iXi ^ 1X2 ^ 1X3 ^ >^1) + {vi V3 ^ V2 ^ Vl) , 


(129) 


with 

f{vi,ix 2 ,i' 3 ) = 640 wg + 48 wg [ 3 i^i -I- dz/f -|- + 41211x2 + 21x11x3 + 41x21x3] 

-t- 4 cuG “t“ 2 z/ 2 “t“ 41 ^ 2^3 F 16 ^ 2^3 F 8r^2^3 F 3r'3 F 4 ix^ {2ix2 F 1x3') 
-|- 2 ^'i { 2 ix 2 F 1x3) {1X2 F 1x21x3 F 21^3) -|- 2 ix^ (51X2 F 5 ^' 2^'3 F 2 z/'|^] 
Fi^l^l (1x2 F 1x3) F 1X1IX2IX3 {1x2 F 1x3) (2 v2 F r'g) 

+ 1 x 1 {1x2 F 2 r'|i /3 -I- Qixlvl + 5 iX 2 ix^ F ^'3) 

+ ixf { 2 v 2 F i^ 3 ) (j^2 F 1X2IX3 F 2 r' 3 ) -|- i>l -|- :/ 2 i ^3 F vl) . ( 130 ) 

It Is a generic feature of the 2PPI approach that, due to the dlagonality of the 
external source J, G Is not frequency dependent; this is neither due to the simplicity 
of the model of the oscillator nor to the truncation. Also, depends on only one 

frequency (instead of three as it would in a 2PI approach). This is a clear advantage in 
the numerical treatment of the problem if indeed only the accessible diagonal quantities 
are desired. 
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Three truncation schemes are investigated in Ref. |46]. The natural one is 


obtained by setting T^"'' to their initial values T^"'' for n > 3. A reduced scheme 

is defined by setting to zero. This reduction allows to gain some analytical 

insights and solve the flow equations analytically [46] 


(n) 


= WG 7 




1 


1 


g/4! ' ^ V’ ^ 2a;G 1 + 5/4!' 

A physically motivated suggestion also leads to an improved scheme [46]: At A = Ai, 
the relation G\ = holds. For X ^ Xi, this relation is inverted to define an effective 
A-dependent frequency characterizing the oscillator w®® = . This is then used 

to calculate according to 


d3).(4) 


= F 


(3).(4) 


,eff 


(132) 


Even without performing explicit calculations, one can state two points: on the one 
hand, the motivation for this ad-hoc replacement is physically comprehensible: it 
corresponds to finding the harmonic oscillator that reproduces the present density at 
the given RG scale A. In a sense, this is similar to the mean-field idea. On the other 
hand, the replacement is unsystematic and carries the danger of double counting of 
diagrams. Whether this is indeed the case, remains to be investigated. Such a study 
is technically not straightforward to conduct and we refrain from pursuing it here. 

Numerical results for the quantum anharmonic oscillator were presented in 
Ref. [15] for finite temperatures T > 0. Our data for the reduced and improved 
DFT-fRG truncation schemes differ from the ones given therej^ As for the other 
fRG schemes, we here only show temperature T — 0 curves. We numerically solved 
the DFT-fRG flow equations using the frequency grid introduced in [Appendix A (the 
reduced scheme has of course an analytical solution). The number of flow equations 
scales ~ mien in the natural as well as in the improved scheme. Taking into account 
the integrations on the right-hand-sides, the effort to calculate the set of differential 
equations in each step of the flow scales ^ n^ien- Results for all three truncation 
schemes are shown in Fig. [^ Note that the main plots show the ground state energy 
and the density as a function of g (and not the relative difference to the exact result 
as above). This is because the reduced and natural DFT-fRG curves are rather far off 
from the exact curve compared to the other RG approaches discussed here. Even for 
this simple problem of the oscillator, these two methods are found to perform much 
worse than simple mean-field theory. In contrast, the improved scheme reproduces 
the exact curve very well. This is discussed in more depth below. 


8. Comparison of the 2PI schemes and conclusion 

In this section, the anharmonic oscillator results obtained from the most successful 
schemes introduced above shall be compared to one another, namely the internally 
closed C-flow, modified D-flow and improved DFT-fRG. Moreover, they shall be set 
into context by comparing them to results of two IPI fRG schemes |1|. These two 

* The first discrepancy is that the free energy curve in the reduced scheme shown in Figure 6 
of Ref. |46| does not agree with the provided analytical result. The second point involves the 
numerical data for the free energy and density calculated in the improved DFT-fRG scheme shown 
in Figures 6 and 7 of Ref. m- The authors of Ref. m informed us that a mistake in their numerical 
implementation is the only origin of the discrepancy to our results. 
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Figure 4. Numerical results for (a) the ground state energy and (b) the density 
as a function of the anharmonicity g. Note that this time the main plots show 
the absolute curves, while the insets show the relative difference to the exact 
curve. We label the schemes by the same names that were used in Ref. |46| for 
better comparability. The natural truncation (DFT-RGS) and improved DFT- 
fRG (DFT-RGI) data are calculated with = 180. For comparison, the exact, 
mean-field (MF) and analytical reduced DFT-fRG (DFT-RG) curves are plotted 
also. 


schemes are based on vertex expansion and can be applied to the fermionic many-body 
problem (cf. Sect.[^ in a straight-forward manner - this is in contrast to fRG schemes 
based on potential expansion. We additionally give brief accounts of the general 
advantages and drawbacks of the different schemes and on the prospects of applying 
2PI fRG to zero- and one-dimensional fermionic quantum many-body problems. 

The first IPI approach to compare with is the regular Matsubara fRG scheme 
that is obtained by truncating the set of flow equations at second order. The second 
is the modification of the first suggested by Katanin [3D]— it corresponds to replacing 
the so called single-scale propagator with a full A-derivative of the full propagator in 
the flow equation for the two-particle vertex. The number of flow equations scales 
^ mfgjj in both IPI fRG schemes. Both schemes use a sharp imaginary frequency 
cutoff as flow parameter. This means that the integrations on the right-hand-sides 
of the flow equations are “canceled” by a (5-function in the regular IPI fRG whereas 
this is not entirely the case in Katanin IPI fRG. Gonsequently, the effort to calculate 
the set of differential equations in each step of the flow scales ~ ^fen regular IPI 
fRG and ~ in Katanin IPI fRG. For the ground state energy and the density 
of the anharmonic oscillator, the original numerical IPI fRG results [I] are used for 
comparison to the 2P(P)I schemes. 

Fig. El shows the numerical results. First of all, all three 2P(P)I schemes 
outperform the IPI ones (except for improved DFT-fRG in the case of the density 
for g > 20). That being said, internally closed C-flow and modified [/-flow yield 
results similar to those of Katanin IPI fRG (but a little better). One should stress, 
however, that modified [/-flow does so with an effort ^ "m-ien each step of the flow 
while the other two methods scale ~ We emphasize that these scaling properties 
also hold for fermionic many-body problems. 

As far as the ground state energy is concerned, improved DFT-fRG produces the 
most accurate results for the shown anharmonicities g. Note however, that similar 
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Figure 5. In this figure, the numerical results of the most successful 2PI 
methods for the ground state energy egs (a) and the density q (b) are compared. 
The relative difference to the exact results is shown. Also, the mean-field (MF) 
curve and the regular IPI fRG (RIPI) and Katanin IPI fRG (KIPI) curves are 
presented. Internally closed C-fiow (ICGF) is calculated at = ^5, modified 

f/-fiow (MUF) and improved DFT-fRG (DFT-RGI) are calculated at m^gj^ = 180. 


to what is observed for the density around g = 20, the improved DFT-fRG energy 
presumably becomes less accurate than the IPI and 2PI ones around g « 60. This 
is strongly linked to the observation that in contrast to the other schemes discussed 
in this section the improved DFT-fRG results for the energy and density deviate 
in a non-monotonic way from the exact one. This behavior might originate from 
the ad hoc replacement ujq —>■ This mean-field-like substitution is at the same 
time at the heart of the success of the improved DFT-fRG scheme. This is plausible 
for the quantum anharmonic oscillator, as mean-field produces rather good results 
for this problem (see above). We expect that the improved DFT-fRG scheme is 
particularly suitable for problems in which the mean-field approximation already 
captures crucial parts of the interaction. In DFT-fRG it is generally impossible 
to access off-diagonal elements of two-point functions. This means that time-non¬ 
local properties of propagators cannot be computed although they are needed for 
the computation of certain observables (e.g. a spectral function). While this is true 
in equilibrium, time-non-local properties play an even more important role in non¬ 
equilibrium situations. Thus, it is not conceivable how to investigate non-equilibrium 
with this scheme. We re-emphasize that improved DFT-fRG might be prone to double 
counting of diagrammatic contributions. Prior to any application to fermionic many- 
body problems this issue should be thoroughly investigated. 

Also the success of the modified Il-How with the Hartree-Fock solution as the 
starting point of the RG procedure heavily relies on the accuracy of mean-field 
theory for the problem at hand. When applied to zero- (quantum dots/impurities) 
or one-dimensional (quantum wires) fermionic many-body problems, modified [/-flow 
faces the problem that the unrestricted Hartree-Fock starting point shows spurious 
symmetry breaking. This raises the interesting question whether the RG procedure can 
restore the symmetry that is broken in the initial conditions. If so, this approximation 
might turn out to be a favorable combination of the aspects of the physics mean-field 
gets right and those which other fRG schemes (such as e.g. IPI fRG [UlSllI]) capture. 
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Figure 6. This figure shows the frequency dependence of the self-energy 
(on a log-linear scale) for several approximate methods at zero temperature and 
for coupling constant g = 10—improved DFT-fRG is missing as the self-energy 
cannot reasonably be computed in that scheme. Internally closed C-flow (ICCF) is 
calculated at = 45 and modified f/-fiow (MUF) is calculated at 


In case the symmetry is not restored by the RG flow, one might resort to restricted 
Hartree-Fock as the starting point or employ plain U-flow for which the effort scales 
^ mfen each step of the flow as well. The internally closed C-flow produces rather 
good results but does not constitute a true advance compared to Katanin IPI fRG as 
it also scales ~ wfen each step of the flow. 

As a last point, the frequency dependence of the self-energy of the anharmonic 
oscillator shall be discussed. Investigating this quantity promises some deeper 
understanding of the Cgsig) and g{g) curves. The frequency dependence for fixed 
interaction g = 10 is shown in Fig. |6|H| All frequency dependent fRG schemes 
underestimate S(^) around = 0; in this regime, the 2PI fRG curves are closer to the 
exact solution. Recall that reproducing the correct behavior for small frequencies is 
more crucial than for large frequencies because T,(v) appears in the calculation of the 
propagator in sum with such that the large ly contribution of S(i^) becomes less 
important. Thus, this observation makes it plausible that the 2PI schemes perform 
better than the IPI ones. 

Acknowledgements 

We are grateful to Sandra Kemler and Jens Braun for numerous discussions. We thank 
the DFG for support via the Research Training Group 1995 ’’Quantum many-body 
methods in condensed matter systems”. 

i As curves for the frequency dependence of the self-energy were not shown in the original work 
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Appendix A. Numerical implementation 
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In our numerical calculations, we set wg = 1- The greatest challenge in implementing 
the zero temperature equations is that the frequency axes are continuous and must 
be discretized for the implementation: v ^ Vn (here, —mien ^ n < mien denotes the 
discretization and not the Matsubara index). We require = —Vn- A cubic spline 
interpolator is used to interpolate function values at v ^ {i^n\ — mien ^ n < mien} 
and to calculate frequency integrals. Almost all right-hand-sides are 

~ 5000 is sufficient to match the accuracy of the differential equation solver. 
Furthermore, it seems advantageous to concentrate the frequencies around zero. 

The following grid definition matches the above criteria (let n > 0): 




nAi 


9i + 


/i 


-1 


/i-i 

nTl ITljen 2 

52/2 


for n < mien,i 

Ah' for mien.i < n < mien ,2 
for mien .2 <n < mien 


(A.l) 


Here, gi,g 2 are externally given parameters. Av = gi/mien,i where integer division 
defines mien.i = mien/3 and mien ,2 = (2mien)/3. /i is determined such that 
^mien 2—92 and /2 is determined such that = 5120. Preliminary calculations 
showed that gi = 2 and 32 = 20 are reasonable choices. Numerical convergence with 
respect to mien was checked carefully. 
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